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We investigate baryon-baryon (BB) interactions in the 3-flavor full QCD simulations with degen- 
erate quark masses for all flavors. The BB potentials in the orbital 5-wave are extracted from the 
Nambu-Bethe-Salpeter wave functions measured on the lattice. We observe strong flavor-spin 
dependences of the BB potentials at short distances. In particular, a strong repulsive core exists in 
the flavor-octet and spin-singlet channel (the 8 S representation), while an attractive core appears 
in the flavor singlet channel (the 1 representation). We discuss the relation of such flavor-spin 
dependence with the Pauli exclusion principle at the quark level. The possible existence of an 
//-dibaryon resonance above the AA threshold is also discussed. 
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1. Introduction 

The generalized nuclear force, which includes not only the nucleon-nucleon (AW) interaction 
but also hyperon-nucleon (YN) and hyperon-hyperon (YY) interactions, has been one of the central 
topics in hadron and nuclear physics. Experimental studies on the ordinary and hyper nuclei as 
well as the observational studies of the neutron stars and supernova explosions are deeply related 
to the physics of the generalized nuclear force 

Concerning the AW interaction, significant number of scattering data have been accumulated 
and are parametrized in terms of the phase shifts or of the phenomenological potentials [||]. On the 
other hand, the properties of the YN and YY interactions are hardly known even today due to the 
lack of high precision YN and YY scattering data. The approximate flavor SU (3) symmetry does 
not necessarily create a bridge between the nuclear force and the hyperon forces because there 
are six independent flavor-channels for the scatterings between octet-baryons. Although there have 
been previous theoretical attempts to fill the gap on the basis of the phenomenological quark models 
[Q] and of the one-boson-exchange models [|J], it is much more desirable to study the generalized 
nuclear force from the first principle QCD. 

Recently a method to extract the interaction potential from lattice QCD simulations has been 
proposed and applied to the AW system [||], and EN and AW systems |§] as well. This paper is 
our exploratory attempt to unravel the flavor-spin structure of the BB interaction on the basis of the 
lattice QCD simulations [Q]. To capture the essential structure, we take the finite but degenerated 
u, d, s-quark masses, so that the system is in the flavor SU (3) limit. In the limit a convenient base 
set exist to look into the structure. Moreover, SU (3) limit result is easy to compare with effective 
models and useful to pin down physical origin of the BB interaction. 

In the next section, we introduce the flavor SU (3) limit and the formalism we employ. In 
section 3, we explain a setup of our numerical simulations. In section 4, we show the resulting BB 
potentials and discuss their implications. The last section is devoted to summary and outlook. 

2. Baryon-baryon potentials in the flavor St/ (3) limit 

In the flavor SU (3) symmetric limit, all baryons and baryonic systems are strictly classified 
into multiplets which serve irreducible representation of the flavor SU (3) group. For example, 
two-octet-baryon systems are classified as 



where "symmetric" and "antisymmetric" stand for the symmetry under the flavor exchange of two 
baryons. The Pauli principle between two baryons imposes, for the system with even (odd) orbital 
angular momentum, 27, 8 V and 1 to be spin singlet (triplet), while 10*, 10 and 8„ to be spin triplet 
(singlet). These six states provide good basis to describe BB interaction. In particular, for 5-wave, 
all off-diagonal interactions vanish and six diagonal one remain. The corresponding potentials are 



808 = 2708,01 ©10*©10©8 a 
v .. ' v .. ' 




(2.1) 



symmetric 



% : V< 27 >(r), V^(r), V«(r), 
3 S, : V( 10 *)(r), V^(r), V^(r) . 



(2.2) 
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Table 1: Summary of lattice parameters and hadron masses. For details about the action and the lattice 



lattice 


P 


a [fm] 


L [fm] 




m ps [MeV] 


m B [MeV] 


A^cfg 


16 3 x 32 


1.83 


0.121(2) 


1.93(3) 


0.13710 
0.13760 


1014(1) 
835(1) 


2026(3) 
1752(3) 


700 
800 



These six potentials contain the essential structure of 5-wave BB interaction. In fact, any potential 
among octet-baryons, including both the diagonal one (B1B2 — > #1^2) and the off-diagonal one 
(B1B2 B3B4), is obtained by suitable combinations of these six potentials. 

According to Refs. [Q], non-local but energy-independent potential can be defined and ex- 
tracted in lattice QCD, from the Nambu-Bethe-Salpeter (NBS) wave function <p(r) and the energy 
E, through the Schrodinger type equation. As the leading order of the derivative expansion of the 
non-local potential, a local potential is obtained by 

where pL is the reduced mass of the system. It is shown that the leading local potential well dominate 
the full potential at low energy region [§]. We use eq. (2.3) and obtain the six potentials in eq. 
(2.2). For 3 Si states, we do not decompose central and tensor potential in this study, and hence the 
resulting potential is the one so called "effective central" potential. 

In the lattice QCD simulations, the NBS wave function $( a \r) for the smallest energy is 
extracted from the four-point function as 

G[ a \t-t ,r) = <0| (BB)W(t,7) W) {a \t ) |0> - 0(«)(?) e -*M>) (2 . 4 ) 



for t — to 3> 1, where (BB) (to) is a wall source operator at time to to create two-baryon states 
in a-plet, while (BB)( a >(t,r) is the sink operator at time t to annihilate two-baryon states. These 
two-baryon operators can be constructed with product of the baryon composite field operator and 
the SU (3) Clebsch-Gordan (CG) coefficients. 



3. Numerical simulations 

We utilize two gauge configuration sets provided by the Japan Lattice Data Grid (JLDG) 
and International Lattice Data Grid (ILDG) 0, which are generated by CP-PACS and JLQCD 
Collaborations JlO|], with the renormalization-group-improved Iwasaki gauge action and the non- 
perturbatively 0(a) -improved Wilson quark action. The lattice parameters and hadron masses are 
summarized in Table [[]. Those hopping parameters K uc j s reside around the physical strange quark 
mass region in 2+ 1 flavor QCD at the same /3 [[[(]] . 

We calculate quark propagators for the spatial wall source at to with the Dirichlet boundary 
condition in the temporal direction at t = to + 16. In order to enhance the signal, we use all 32 time 
slices on each configuration as to, and take the average over forward and backward propagations in 
time. We estimate the statistical errors by the jackknife method. All these numerical computations 
have been done at the KEK supercomputer system, Blue Gene/L and SRI 1000. 
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Figure 1: Lattice QCD extracted potentials V^ 27 ' in the left and y( 10 *' in the right. 



4. Results and their implications 

4.1 BB potentials in flavor basis 

Figure |l| to |5| show the obtained six BB potentials. Red bars (green crosses) correspond to 
the pion mass of 1014 MeV (835 MeV). Since the effective mass of the four-point function in 
each channel shows a plateau at t — to > 10, we use data at t — to = 10 exclusively to extract both 
NBS wave functions and BB potentials throughout this paper. We have estimated the energies 



contribution (the second term on the right-hand side of Eq. ( J2.3[ )) under the condition that Eq. (2.3) 
approaches zero at large r (the method called "from V" in Ref. [||]). We find £"( 27 ) ~ —5 MeV, 
£(*•) ~ 25 MeV, EW ~ _30 MeV, £( 10 *) ~ -10 MeV, E^ ~ MeV and E^ ~ -15 MeV. 
Since the energy contribution is small compared to the Laplacian term (the first term on the right- 
hand side of Eq. (2.3)) at a short distance, we have not attempted to extract a precise value of E^ 



in this study. 

Figure |l] shows potentials W 27 ) and V^- 10 '. Because a symmetric AW belongs to 27 and anti- 
symmetric AW belongs to 10* in the St/ (3) limit, V (27 ) and V (1(r) can be regarded as the St/ (3) 
limit of the nuclear force potential in l So and 3 Si, respectively. Both potential have a repulsive 
core at a short distance and an attractive pocket at around 0.6 fm. These qualitative features agree 
with our previous results found in the AW system in quenched approximation with lighter quark 
mass[||]. This gives good consistency check. Furthermore, these feature qualitatively agree with 
well known property of the phenomenological AW potentials constructed form data. This shows 
efficiency of the present method. 

Figure |2| shows potentials and V^ 8a \ Since anti-symmetric AX(/=3/2) belongs to 10 in 
the SU (3) limit, one can regard W 10 ) as the SU (3) limit of 3 Si AX(/=3/2) potential, for example. 
Similarly, is the St/ (3) limit of 3 Si AtE(/=0) potential. As you see, has a stronger 
repulsive core and a much shallower attractive pocket than yf 2710 ). On the other hand, y( 8 °) has 
a very weak repulsive core and a deeper attractive pocket compared to y( 27 > 10 ). This is exactly a 
part of the essential structure of BB interaction we looked for. It is interesting to point out that the 
above features were foreseen in a simple quark model, where effect of quark Pauli blocking and 
one-gluon-exchange interaction are evaluated systematically using the SU (6) quark wave function 
of baryon^]. We come back this point later. 
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Figure 2: Lattice QCD extracted potentials V^ 10 ) in the left and V ' 8 "' in the right. 
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Figure 3: Lattice QCD extracted potentials V^ 8 - 1 ' in the left and VW in the right. 



Figure |3] shows potentials and V"W. These potentials cannot be interpreted as an SU(3) 
limit of some potential between two baryons, but are always superpositions of them. First, we see 
that has a very strong repulsive core. It is strongest among all the channels. On the other 
hand, shows no repulsion at a short distance but attraction instead. It is in contrast to other 
cases. Now, the essential flavor-spin structure of the S-wave BB interaction is revealed completely. 

As the case of V^ 10 ' and V^ a \ the above features of W 8 " and V^ 1 ' were also foreseen in 
simple quark models^, [jl]]. In particular, in quark models becomes strongly repulsive at a 
short distance since the six quarks cannot occupy the same orbital state due to quark Pauli blocking. 
On the other hand, in the 1 channel six quarks do not suffer from the quark Pauli blocking at all, 
and the potential can become attractive due to short-range gluon exchange. In total, simple quark 
models predicted the essential feature of BB interaction correctly. This fact suggests that the quark 
Pauli blocking plays an essential role in the repulsion in BB system. 

4.2 BB potentials in baryon basis 

With explicit flavor SU (3) breaking, the baryon basis such as AW and AN are more suitable 
to describe BB interactions than the flavor basis. Generally, interaction potential takes matrix form 
in coupled channels in the baryon basis. In the SU (3) limit, the baryon-base potentials can be 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



r [fm] t [fm] 

Figure 4: BB potentials in baryon basis for S=— 2, 7=0, 1 Sq sector. Three diagonal (off-diagonal) potentials 
are shown in the left (right) panel. Phases of off-diagonal ones in the right panel are arranged in a zoom-out 
plot. Their true signs are shown in the insert. 

obtained by a unitary rotation of flavor-base potentials as V, ; (r) = £ a Ui a V^ a '(r) U a , where U is a 
unitary matrix made of CG coefficients, which rotates the flavor basis to the baryon basis. 

In Figure 0, we show the potentials for ^o, S=— 2, 1=0 sector as characteristic examples. To 
obtain Vj/(r), we have used a simple analytic function with five parameters fitted to data of V^ a \ 
The left panel of Fig. ^ shows the diagonal part of the potentials. We see that AA and ,/VS(/=0) 
channels are attractive while LL(I=0) channel is repulsive. The attraction in the 1 channel is most 
reflected in the jVS(/=0) channel, while the strong repulsion in the 8 S channel is most reflected in 
the ££(7=0), owing to its largest CG coefficient among the three channels. Nevertheless, all three 
diagonal potentials have a repulsive core originating from the 8 V component. The right panel of 
Fig. f| shows the off-diagonal parts of the potential matrix, which are comparable in magnitude to 
the diagonal ones. Since the off-diagonal parts are not negligible in this sector, full coupled channel 
analysis is necessary to study the observables. 

Similarly, we have obtained potentials for all other sectors with different strangeness, iso-spin 
and spin, and find their features^]. Since those features may remain in the flavor-»S£/(3)-broken 
physical world essentially, we expect our result provide useful hints on the behavior of hyperons 
(A, £ and S) in hypernuclei and in neutron stars 

5. Summary and outlook 

We have performed 3-flavor full QCD simulations to study the general features of the BB 
interaction in the flavor SU (3) limit. From the NBS wave function measured on the lattice, we 
have extracted all six independent potentials in the S-wave at the leading order of the derivative 
expansion. We have found strong flavor-spin dependence of the BB interactions. It turns out there 
exist variety of BB interaction in u, d, s-quark three-flavor world, in contrast to u, d-quark two- 
flavor world. In particular, V^ 8 ^ has a very strong repulsive core at a short distance, while is 
attractive at all distances. These features are consistent with the Pauli blocking effect among quarks 
previously studied in phenomenological quark models. This indicate that the quark Pauli blocking 
is essential for the repulsion in BB interactions at a short distance. 
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In the flavor singlet BB channel, a bound state so called //-dibaryon is expected in phenomeno- 



logical quark models Qll[]. The present lattice QCD data are insufficient to derive a definite con- 
clusion on the //-dibaryon, because of a single and small lattice volume. Now we are performing 
lattice QCD simulations at flavor-S£/ (3)-point in larger volume. Nevertheless, using the potential 
V^\r) and solving the Schrodinger equation, a shallow bound state is found, suggesting a pos- 
sibility of a bound //-dibaryon with baryon components AA : EE : EN = — vT : v3 : y4, in the 
SU (3) limit world with heavy quark. In the real world, flavor SU (3) breaking is considerable and 
the thresholds of the S=— 2, 1=0 sector are located in order AA < EN < EE. In our test study 
of solving the Schrodinger equation with small baryon mass differences from a 2+1 flavor lattice 
QCD simulation, a resonance state is found at an energy between AA and EN threshold. If such a 
resonance exists in nature, it may explain the enhancement just above the AA threshold recently re- 



ported in the KEK experiment [12Q. However, further investigations in both theory and experiment 



are necessary to derive a definite conclusion. 
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